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ABSTRACT 

The stratification near the base of the Sun's convective envelope is governed by processes of 
convective overshooting and element diffusion, and the region is widely believed to play a key 
role in the solar dynamo. The stratification in that region gives rise to a characteristic signal 
in the frequencies of solar p modes, which has been used to determine the depth of the solar 
convection zone and to investigate the extent of convective overshoot. Previous helioseismic 
investigations have shown that the Sun's spherically symmetric stratification in this region is 
smoother than that in a standard solar model without overshooting, and have ruled out simple 
models incorporating overshooting, which extend the region of adiabatic stratification and 
have a more-or-less abrupt transition to subadiabatic stratification at the edge of the overshoot 
region. In this paper we consider physically motivated models which have a smooth transition 
in stratification bridging the region from the lower convection zone to the radiative interior 
beneath. We find that such a model is in better agreement with the helioseismic data than a 
standard solar model. 
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1 INTRODUCTION 

An understanding of the overshoot region at the bottom of the Sun's 
convective envelope is important for a number of reasons. The over- 
shoot region approximately coincides with the solar tachocline, a 
region of rotational shear which is generally believed to play a key 
role in the solar dynamo: overshooting is likely to be important for 
helping to store the magnetic flux below the convection zone during 
the solar cycle. Bulk motion in the overshoot region also affects the 
thermal stratification and it may contribute to significant mixing of 
chemical elements, for example transporting fragile elements such 
as lithium to hotter regions where they are destroyed more easily 
than in the convection zone. More generally, convective overshoot 
in stars (particularly those with convective cores) is likely to be an 
important and as yet imperfectly understood process affecting age 
estimates of stars, and so improved constraints on theories of over- 
shooting obtained from a study of how overshooting works in the 
solar case can be important for understanding stars more widely. 
Helioseismology provides a means of probing directly the 
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conditions inside the Sun, because the frequencies of resonant 
modes set up by acoustic waves propagating in the solar interior 
depend in particular on the adiabatic sound speed c which is given 
by 

P m^n 

Here Fi is the logarithmic derivative of pressure p with respect to 
density p at constant specific entropy, T is temperature, ji is the 
mean molecular weight, is Boltzmann's constant and ot„ is the 
atomic mass unit. Hence the sound-speed gradient with respect to 
depth depends on the temperature gradient, which itself depends 
on the mechanism by which heat is transported. The transition be- 
tween fully radiative heat transport beneath the convection zone 
and convective heat transport within the convection zone is mani- 
fest in the temperature gradient and hence too in the sound-speed 
gradient. If the transition in sound-speed gradient takes place over 
a distance that is small compared with the vertical wavelength of 
the acoustic waves near the base of the convection zone, then the 
transition appears to the waves to be more-or-less sharp and this 
gives rise to an oscillatory signal in the mode frequencies oj: the 
form of the signal gives information about the location and nature 
of t his "acoustic glitch" ItG ough 2002a). 

iMonteiro et alj ( Il994h represented the efi"ect of the base of the 
convection zone in terms of an additive contribution 6(i>o to the 
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frequencies, relative to those of a corresponding model in which 
the transition had be en smoothed out. For low-degree modes (see 
iMonteiro et aDl2000h a break in the first derivative of the sound 
speed gives rise to a contribution of the form 



6a)„ = A(w) cos (2(Df(i + 20o) , 



(2) 



while a break in the second derivative gives rise to a similar signal 
but with a sine term instead of a cosine, and with a different fre- 
quency dependence of A(w). Here fd is essentially the value of t at 
the location of the acoustic glitch, where 



(3) 



is acoustic depth beneath the surface, r being the corresponding dis- 
tance to the centre and R the surface radius of the Sun. Also, tpg is a 
phase introduced by the reflection of the mode at the turning points, 
depending in particular on the near-surface structure. The function 
A(ci)) is an amplitude which depends on the sharpness and nature of 
the convection-zone base: the smoother the transition, the smaller 
in general will be the amplitude. However, if moderate-degree data 
are used, as in the case of Sun where we have accurate data for 
modes whose degree is above 3, the above expr ession needs to in- 
clude additional terms, both in the amplitude (Monteiro et al. 1994") 
and in the argument of the signal (Christensen-Dalsgaard et al., 
Il995h . to account for the first-order effect of the mode degree on 
the signal. 

Overshoot at the base of the solar convection has traditi onally 
been modelled using non-local mixing-length theory (e.g.. IZahg 
Il99lh . Such models mostly predict an overshoot region that is 
nearly adiabatically stratified; in terms of the logarithmic tempera- 
ture gradient V = din T/dln p one finds that 6 = V - Vad 10"*, 

where Vad = (31n T/cJlnp), is the adiabatic temperature gradient, 
s being specific entropy. The depth of overshoot region is typi- 
cally between Q.lHp and QAHp, where Hp is the pressure scale 
height at the base of the convection zone, with a very steep tran- 
sition towards the radiative temperature gradient. The results have 
seemed rather robust, since the above behaviour is found in models 
inc orporating quite d i fferen t large-scale flow structures: for exam- 
ple |yan^aiiegooije3 ^^13) assumed overturning convective rolls 
whilst Schmitt et al. r( ll984h explicitly modelled downward plumes 
in the overshoot region. We note, however, that other treatments 
of the overshoot region have suggested a much smoother transition 
to the radiative gradient (e.g., Xiong & Deng||200ll : lDeng & Xiond 



one tenth of a pre ssure s cale-height. A similar limit wa s found by 
iBasu et al.l h994h . while lRoxburgh & Vorontsovl ( 1 19941) obtained a 
somewhat weaker limit. iBasu & Antial ( Il994l) mated that the com- 
position gradient caused by the inclusion of helium settling pro- 
duced a sharper transition in the sound-speed gradient, even in 
models without overshoot, and hence a larger oscillatory signal. 
From these analyses it would appear that the transition in sound- 
speed gradient at the base of the solar convection zone is if any- 
thing smoother than in the non-overshoot models. A caveat is that 
what we purport to measure in the above studies is the spherically 
symmetric component of the structure: departures from sphericity, 
such as a latitudinal dependence to the shape of the base of the 
convection zone, could make the transition appear smoother than 
it is locally. The helioseismic evidence, however, is that the loca- 
tion of the base of the con vection zone is independent of latitude 
jMonteiro & ThompsorJl998, : Basu & Antia 2001.) . Changes of the 
base of the convection zone on time scales shorter than the obser- 
vation interval could also have a similar effect by introducing a 
time- averaging effect on the mode frequencies that would mimic 
a smoother transition. The importance of this effect is difficult to 
estimate, however, and depends strongly on the 3D nature of con- 
vection at the base of the envelope. 

Overshoot has been addressed in the last two decades by a va- 
riety of 2D and 3D numeri c al simulations CRoxburgh & SimmonsI 
19931 : iHurlburt et al.l Il994 ISingh et all ll995L Il998l: ISaikia et al.l 



20od : iBrummell et al.l |2002| : iRogers & Glatzmaierl l2005bllal) 
Whilst the non-local mixing-length models have clearly predicted 
an adiabatic overshoot region of a sizeable fraction of a pressure 
scale height and a rather sharp transition to the radiative zone be- 
neath, the numerical simulati ons show a greater vari ety of possi- 
ble behaviours. The work by IBrummell et al.l ( |2002|) is currently 
one of the best resolved and most turbulent investigations: it shows 
strongly subadiabatic overshoot with very smooth transition to- 
wards the radiative temperature gradient. Most of the earlier 2D 
and 3D simulation were more in the laminar regime and found, 
depending on their parameters (mainly the stiffness of the subadia- 
batic layer), both nearly adiabatic overshoot and extended subadia- 
batic overshoot. 

iRempell ( l2004h tried to understand the cause of the discrepan- 
cies in different overshoot treatments, which have been attributed 
either to over-simplifications in the mixing-length approach and re- 
lated models or to the fact that numerical simulations are not in 
the correct par ameter range. Using a semi-analytic plume model 
lRempeil ( l2004l) showed that the main differences result from differ- 
ent values of the energy flux used in these models (more exactly the 
energy flux divided by the filling factor of downflows, expressed by 
a dimensionless number O = F/(f p ^jplg), with the energy flux F, 

. . . .. the downflow filling factor /, pressure p and density g at the base 

1994 IChristensen-Dalsgaard et al.lll995 ) in order to compare the of the convection zone). The value of O ranges from 10 (mixing- 
length models) to up to 10"^ (numerical simulations). The influence 
of the energy flux (lower energy flux leading to more adiabatic 
oversh oot) was already indicated i n the work of Brurnmell et al.l 
( I2OO2I) and has been confirmed by iKapvla et alj hoOTi) by vary- 
ing the energy flux over two orders of magnitude (with the caveat 
that lowering the energy flux in a numerical simulations typically 
decreases al so the degree o f turbulence). The main conclusion of 
the work of lRempell ( |2004|) is that non-local mixing-length mod- 
els and current numerical simulations mark two extremes in terms 
of the parameter O and the Sun might have a solution somewhere 
in between, however most likely closer to the mixing-length regime 
unless the downflow filling factor at the base of the convection zone 
is tiny. This leads to the possibility that solar overshoot has a strat- 



I2OO8I : iBaturin & Mironovair2010l) 

The rather abrupt transition in the temperature gradient pre- 
dicted by the non-local mixing-length models has been param- 
eterized and incorporated into solar models ( Basu et aU 19941: 
Basu & Antialll994l:lMonteiro et aPl 1994 [Roxburgh & Vorontsovl 



predicted acoustic-mode frequencies with those observed on the 
Sun. Models with and without overshooting all have an oscilla- 
tory signal in frequencies coming from the base of the adiabatically 
stratified region, but the amplitude of the signal in the frequencies 
in the overshoot models is greater than in models without over- 
shooting. 

When the observed and model frequencies are compared, it 
is found that the amplitude of 6lo^ in the Sun is comparable with 
or smaller than that in models without overshooting, implying 
that the amount of overshooting of the kin d predict e d by these 
mixing-length models is ve r y sma ll. .Monteiro et al. I ( 1 1994 and 
IChristensen-Dalsgaard et al] i ll 9951) concluded that the extent of 
any such overshoot at the base of the convection zone was less than 
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ification close to adiabatic, but with a transition towards the radia- 
tive gradient that is smoother than found in non-local mixing-length 
models. 

Accordingly, we have created a family of semi-analytic mod- 
els that display a range of possible characteristics for the convection 
zone, which may give better agreement with the helioseismic ob- 
servations than the old non-local mixing-length models. Moreover, 
helioseismology may determine which model provides the best fit 
to the actual stratification of the Sun; this may, we hope, teach us 
something about the physics of convection and convective over- 
shooting in stellar interiors. 

The goal of the present paper is to discuss these models and 
compare their seismic properties with those observed. Section 2 
presents the characteristics of the semi-analytic overshoot model 
and its parameterization in our stellar structure calculations. In Sec- 
tion 3 we describe the seismological method we use to analyse the 
oscillation frequencies to establish the characteristics of the transi- 
tion in structure near the base of the convective envelope, and we 
apply it to the analysis of solar data in Section 4. Section 5 describes 
the solar models used in this study and presents the analysis of the 
frequencies computed from those models. In Section 6 we synthe- 
size the observational and model results and discuss our inferences, 
and we present our conclusions in Section 7. 



0.45 




Figure 1. The temperature gradients V (solid curve) and V^d (dotted curve) 
in Model C2 with overshoot characterized by f'cz - ''i = 0.01, r^- rf = 0.02 
and p = 0.5. The thin vertical dashed lines show the location of (from left to 
ri ght) Tf, Ti and r^x- For compar i son, t he dashed curve shows V in Model S 
of^ lChristensen-Dalsgaard et al.l i ll 9961) . with no overshoot. The horizontal 
bar shows the range, around rt, over which the argument ijj (cf. Eq.[8), at a 
frequency cojln = 2500 /jHz, changes by n and hence provides a measure 
of the vertical wavelength of the acoustic waves. 



2 A SIMPLE PARAMETERIZATION FOR OVERSHOOT 

As summarized in the introduction a variety of ov ershoot models 
including simplified plume mod els (Rempel 2 004|). non-local con - 
vection models (e.g., IXiong"& Dena,2001.: ,Deng & XiondbOOSl) , 
and 3D simulations (e.g.. lBrummell et alj|2002h predict overshoot 
profiles that are substantially smoother than those obtained through 
non-local mixing-length models. While the physical reasons for the 
smoother transitions dilfer among these models, the resulting pro- 
files show a large degree of similarity. For the purpose of the helio- 
seismic investigation in this paper it is not essential to use exactly 
one of these models. Rather what is important is to realize that over- 
shoot profiles smoother than those obtained by standard mixing- 
length theory are possible. Using the models o f iRempell (2004) as 
guidance we find the following properties to be rather general for 
all profiles: 

(i) The overshoot profile matches smoothly with the radiative 
gradient beneath. 

(ii) In the overshoot region V is in between Viad and Vaj. 

(iii) The lower part of the convection zone is weakly subadia- 
batic. 

In this investigation we use a parameterization which produces 
overshoot profiles possessing the above-mentioned properties. Af- 
ter finding the overshoot profile that is most consistent with helio- 
seismic data we shall return to overshoot models and discuss po- 
tential implications. 

The idea is to provide an analytical match to the behaviour ob- 
tained in the numerical simulations by Rempel (2004), constrained 
such that the temperature gradient and its first derivative are every- 
where continuous. This is done by representing the actual temper- 
ature gradient as a function of r thus: 



V = Vad - T„v 



for r ^ Tf . 



(4) 



7? + exp(2<r) 

where ( = {r - rt)/d. For r < rf, the temperature gradient is ra- 



diative, V 
that 



V = V 



The constants Toys and d are determined such 



at 



(5) 



where the dash indicates differentiation with respect to r. Thus the 
formulation is characterized by r,, which determines the overall lo- 
cation of the transition from the adiabatic to radiative temperature 
gradient, and rf, which is the radius at the bottom of the overshoot 
region, r,-rf essentially controlling its width. We assume the full 
overshoot region, down to rf, to be chemically fully mixed. The 
parameter p provides additional flexibility to the location of the 
overshoot region, relative to the base, r = r^-^, of the convectively 
unstable region, defined by V^d = Vad. 

Given n, rf and j8, Tovs and d are determined from Eqs (|5), 
which yield 



^v. = - [Vad(rf) - V,-ad(rf)] [/3 + exp(2ff)] 
where ff = (rf - ri)/d, and, neglecting V' ,, 



(6) 



1 2exp(2Z'f) r n 

on using Eq. l|6j. Equation ^ can be solved for d and Eq. ^ then 
yields Toy,,. 

An example of the resulting V is illustrated in Fig. [T] com- 
pared with the reference Model S with no overshoot. This clearly 
illustrates the region of subadiabatic stratification in the lower part 
of the convection zone, and the smooth match to the radiative gra- 
dient at the bottom of the overshoot region, in the model with over- 
shoot. For comparison, the horizontal bar shows the wavelength of 
the squared vertical displacement eigenfunction f (cf . Eq. [8]l at a 
reference frequency a)/2n = 2500 yuHz; it is evident that in both 
models the dominant transition takes place over a distance substan- 
tially smaller than the wavelength. 
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3 METHOD OF SEISMIC ANALYSIS 

Follow ing iMonteiro et alj ( 1 19941) and lChristensen-Dalsgaard et alj 
(Il995h . we consider the frequencies of the modes to be a sum of 
two components: a smooth component u)q, and an oscillatory com- 
ponent ScDp coming from the acoustic glitch caused by the base of 
the convection zone and any overshoot. 

The oscillatory signal arises because of the variation of the 
phase of the mode eigenfunctions at the location of the sharp fea- 
ture, as a function of mode frequency and to a lesser extent of mode 
degree. Asymptotically, the radial displacement eigenfunction is 
given by 

oc (pcY'l^--' cos[i/.(r)] , 



where 



(8) 



(9) 



and is a phase function (it is principally a function of frequency) 
which depends on conditions near the surface of the Sun: we dis- 
cuss if) further below , and in de tail in the Appendix. 

As shown bv lChristensen-Dalsgaard et al.l ( 119951) , the acoustic 
glitch formed by the base of the convection zone and any region of 
overshooting gives rise to an oscillatory signal, with argument 



A(co, I) = IcDTi - 7d + 200 , 

which can be expressed as 
1-2A 



SiL), 



fli |— j sin[A(w, /)] 
1-2A /&\ 

— — — ^2 - cos[A(a;, I)] . 
(1-A)- \cdI 



" (l-A)5/2 
1-2A 

-I- 



(10) 



(11) 



The second term is the form expected from a discontinuity in the 
first derivative of sound speed, and the first term from a discontinu- 
ity in the second derivative. Here 



/(/+1) (cb\- 



A = ^rr-^ - Ad 

/(/+1) 



(12) 



(Aj is defined below) and / and ci) are reference values of / and to. 
For the present investigation we chose the reference values /=20 
andiIi/27r=2500pHz. 

To obtain Eq. jl lb an expansion in A has been used to derive 
the expression for the amplitude: moreover, it has been assumed 
that the acoustic glitch is adequately represented by a single loca- 
tion to perform the expansion in frequency. 

The phase function 0, due to the reflection of the modes at 
the inner and outer turning points, is assumed to be represented 
asymptotically jMonteiro et aljl994h by 



0((jj, /) ^ -I- a^dj + a. 



1(1+1) 
2co 



(13) 



where (a^,ay) are unknown expansion coefficients. The remaining 
parameters are related to the internal structure of the star through 
the following relations: 



r 



dr 



7d 



r 



dr -I- 



/(/+ 1) (c 



(14) 



while Oi and fl2 depend on the structure at the transition (see 
IChristensen-Dalsgaard et al]|l995L for the details). Note that the 
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Figure 2. Difference in zero-frequency phase between the observations and 
Model S , obtained from differential asymptotic analysis (see Appendix A), 
plotted against the mean date of the obsei'vations. 



measurable quantity fd is not identical to tj, the acoustic depth of 
the acoustic glitch, because it includes a contribution from the 
dependence of (/> on frequency. Similarly, yd contains a contribution 
from the /-dependent part of 0. 

In the present work we are mainly interested in exploring the 
nature of the profile of the transition near the base of the convec- 
tion zone. Thus we wish to measure the difi'erent contributions to 
the signal coming from the discontinuities in the first and second 
derivatives of the sound speed. 

In previous works we have found that if (po is also fitted as 
a free parameter, there is a strong correlation between the result- 
ing values of fd and (f>o. This probably results from the fairly nar- 
row range of frequencies included and the weak dependence on 
mode degree for most of the points, which means that a change 
in the average of tufi can be compensated by a change in 0o and 
Ad. To overcome this difficulty the two terms in expression ( II It 
were previously combined in a single cosine function with a single 
amplitude that depends on frequency and mode degree. However, 
this option weakens our capacity to study the actual behaviour of 
ai and a2 for difi'erent profiles of the transi tion. Thus, in a depar- 
ture from what we have done previousl y in iMonteiro et al.l (Il994h 
and Ichristensen-Dalsgaard et aP d 19951) . here we use directly the 
expression dllt in the fit to the data to determine the four parame- 
ters (fli, (32, td, 7d)- The two components of the signal are fitted in 
the frequencies of modes with degrees 5 < Z < 20 and cyclic fre- 
quencies 1800 < LL)/2n < 4000 yuHz using an iterative least-squares 
procedure. The numerical approac h is the same as described in Ap- 
pendix C of lMonteiro et alj ( 119941) . 

As noted above, if 0o is fitted as a free parameter it may be 
highly correlated with fd. However, as discussed in Appendix A, 
00 can be determined independently from analysis of the oscilla- 
tion frequencies and hence can be regarded as known in the fit to 
Eq. nil . We therefore adopt this procedure. There is still the diffi- 
culty of having Ad as free parameter of the fitting, due to the con- 
nection with 7d- Considering that Ad is a slowly varying function 
of T in all models, we have chosen to adopt a reference value from 
the standard solar model corresponding to Ad = 0.35. 

In order to compare the amplitude of the signal for different 
models and the solar data it is convenient to use a reference value 
of the amplitude at fixed frequency O) and at / = defined by 



1/2 



(15) 
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The principal quantities characterizing the signal due to the 
transition near the base of the convection zone are the period, mea- 
sured by fd, and the amplitude (given by a\ and 02). The secondary 
quantities and are necessary to account for the dependence 
on mode degree of t he amplitude and the period of the signal (see 
Fig. 3 of lMonteiro et a l. 1994) as data up to / = 20 is used. The ad- 
vantage of including the higher-degree data is that in general their 
frequencies are more precisely determined, and the much higher 
number of mode frequencies being used renders the fitting much 
more reliable. 

Our least-squares fitting procedure weights all data equally, 
i.e., it does not take the quoted observational uncertainties into ac- 
count. The reason for this choice is that we find it gives better deter- 
minations of the two amplitude parameters ai and 02: otherwise the 
largely systematic variation in mode frequency uncertainties with 
frequency causes these parameters to be less well constrained in 
the fitting. In practice it would correspond to putting all the weight 
on the low-frequency range of the spectrum. This would mainly 
lead to the determination of one of the coefficients, rendering the 
other unnecessary in the fit. The solution would then be mostly de- 
teiTnined by the initial guess. Such a frequency-dependent bias of 
the fitting can be removed if the fitting uses equal weights. 



4 ANALYSIS OF SOLAR DATA 

The solar data used in this work have been obtained from 
MDI/SOHO observations over a total period of over 14 years. To 
determine the signal we use set s of frequenc iej] obtained for 72 
days of continuous observations ( ISchoulll999l) . Sets of frequencies 
for 1-year periods are also considered. These have been obtained by 
combining 5 consecutive sets of 72-d frequencies to produce an av- 
erage for each 1-year set. Finally we considered frequencies from a 
coherently analysed set of 1-yr data from the beginning of the MDI 
observations. The results from this set were very similar to those of 
the corresponding 1-yr average of the 72-d sets and hence will not 
be discussed separately in the following. 

We use modes having degree 5 < / < 20 and cyclic frequency 
between 1900 and 4000 yuHz. Because the least-squares fitting ig- 
nores the observational uncertainty on the frequencies, only modes 
with a quoted observational uncertainty below O.lyuHz are used. 
The phase constant <po was estimated as discussed in Appendix A, 
from a differential asymptotic analysis of each frequency set rela- 
tive to the frequencies of Model S that was adjusted to match ap- 
proximately the solar frequencies. The resulting phase offset S4>o, 
as a function of the date of the observations, is illustrated in Fig.[2] 
In the analysis this offset is added to the phase (fiQ = 0.5485 ob- 
tained for Model S to obtain the appropriate phase constant for the 
particular observational dataset. The variation of 6(f>o with phase in 
the solar cycle is discussed in Section |A3l 

To illustrate the quality of the fitting. Fig. [3] shows the sig- 
nal ScOp for the case of frequencies obtained from a 1-year time- 
series of observations from MDI. In order to show that our pro- 
cedure has cleanly separated the smooth and oscillatory compo- 
nents of the frequencies, the lower panel shows the second dif- 
ferences of the "smooth" component of the frequencies, obtained 
after removing from the solar frequencies the fitted value of the 
signal at each frequency and mode degree. This illustrates that the 



' The frequencies are obtainable online courtesy of J. Schou at 
|http : / /quake . Stanford . edu/~schou/anavw72z/| 
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Figure 3. Top panel: Signal for a set of frequencies from MDI (1-year data), 
shown as filled circles for 5 < ^ < 14 and open circles for 15 < ^ < 20, 
with 1 o" error bars. Also shown are the fitted values at each point from 
the expression in Eq. i ll ft , using the parameters obtained by fitting the fre- 
quencies; to guide the eye, all points have been ordered by the reduced fre- 
quency used as abscissa and connected by a fine. (The somewhat in'cgular 
behaviour aiises from the dependence of the amplitude on a) and /.) Bottom 
panel: Second difterences calculated for the smooth component cjq of the 
frequencies, obtained by subtracting from the frequencies the values rep- 
resented by the line in the top panel. The lines connect points of the same 
mode degree (symbols as above). 



signal coming from the base of the convection zone is no longer 
visible, while the signal due to the helium ionization zone (e.g. 
'Mon teiro & ThomDsonll200^ can be clearly seen in the second dif- 
ferences. 

The impact of the observational uncertainties on fj, A2.5 and 
the other measured parameters has been estimated in two ways. 
Firstly we compute the standard deviations of the measured values 
for a given parameter amongst all the 72-d observational datasets. 
Secondly, we make 500 Monte Carlo simulations of frequency 
datasets, using the estimated errors for a set of 72-d observational 
frequencies, fitting those artificial datasets in the same way as we 
fit the real data, and calculate standard deviations for the resulting 
estimated parameters. The estimate based on all the observed fre- 
quency sets will be affected by any temporal variations in the solar 
values, since the data are obtained at different epochs. The sec- 
ond estimate is purely an indicator of how the uncertainties in the 
frequencies impact the parameters measured from the signal. Both 
sets of uncertainties are shown in Table [T| together with the aver- 
age values for the Sun. It is evident that the standard deviations in 
the parameters inferred from the observations are largely consistent 
with the result of the Monte Carlo simulation. Thus the scatter in 
the inferred values for A2.5 for the observations is essentially con- 
sistent with the assumed error in the observed frequencies. 

Figure |4] shows the infeiTed amplitude A2 5 and period fd for 
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Figure 4. Values of the amplitude /I2.5 (top panel) and acoustic depth fj 
(bottom panel) of the signal for sets of 72 days of MDI frequencies (full cir- 
cles) and for sets of 1-year MDI frequencies (open circles), plotted against 
the mean date of the observations. The horizontal error bar represents the 
period of the time-series used for the calculation of the frequencies. The 
horizontal dashed line corresponds to the mean, and the horizontal dotted 
lines indicate ±3 <t standard deviations obtained by Monte Carlo simula- 
tions based on 72-d uncertainties of the solar frequencies. 



Table 1. Seismic parameters (averages) resulting from the fit of the signal 
in Eq. (TT) to the solar frequencies. Data sets for 72 days have been consid- 
ered, as well as a coherently analyzed 360-d data set together with 1-year 
averages of 5 consecutive sets of 72-d frequencies. The values CTobs are ob- 
tained as the 1 o" distribution of the values found for each group of sets of 
solar frequencies (72-d and 360-d) while <tmc sie the lo" of the Monte Carlo 
simulation (500 realizations) for a single set of frequencies using the quoted 
observational uncertainties in the frequencies for one of the 72-d sets and 
for one of the 1-year sets. 
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Figure 5. The filled circles are the values of A2.5 (upper panel) and fa (lower 
panel), versus fj, for the solar data (72-d frequency sets). Also shown, as 
stars, are the average values of these values together with the error bars 
corresponding to a 1 tr uncertainty estimated from Monte Carlo simulations 
for the uncertainties of a 72-d set of solar frequencies (see Table[T). 



the solar data, as functions of epoch. The dashed lines indicate the 
average values over the whole period, and the dotted lines indicate 
the 3 (T error interval obtained from Monte Carlo simulations based 
on the eiTors in the 72-d datasets. When 1-year data are used (see 
Table[Tll, there is some indication that the variation in the Sun in 
and is slightly above the noise of the data. This could point to- 
wards a variation in time of the stratification of the overshoot region 
or the near-surface layers; there is indeed perhaps some hint of an 
orderly systematic variation of the 1-yr results in Fig. 2] (although 
not at an 1 1-yr period), but more accurate data would be necessary 
to confirm this possibility. 

As a further indication of the properties of the parameters in- 
ferred from the 72-d datasets, Fig. |5] shows A2.5 and against td. 
While A2.5 and f j show no obvious correlation, there is a clear cor- 
relation between and f^. In fact it may be seen from Eq. dl It that 
both f d and fi affect the phase of the signal, in such a way that an 
increase in can be compensated by an increase in yd, as observed 
in the lower panel of Fig. |5] This correlation is also confirmed by 
Monte Carlo simulations carried out on the basis of the model fre- 
quencies (see below). For other pairs of parameters no significant 
correlation was found. 



5 ANALYSIS OF SOLAR MODELS 

The signal for the models has been fitted using the same range of 
mode degree and frequency as done for the Sun (see above). In this 
way we have ensured that the fitting takes place under the same 
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conditions as for the solar data, both for error free data and for the 
Monte Carlo simulations used to estimate the impact of noise. 

5.1 Models emulating overshoot 

Apart from the treatment of the overshoot region, the models 
co nsidered here were compute d in th e same way as Model S 
of IChristensen-Dalsgaard et aL The computation used 

the OPAL equation of state ( Rogers et al. 1996) and opacity 
( Ilglesiasetal.l ll992h , the Bahcall & Pinsonn eaulTl 19951) nuclear 
reaction rates, and the lMichaud & Proffitli H i 9931) treatment of dif- 
fusion and settling of helium and heavy elements, the latter be- 
ing treated as fully ionized oxygen. All models were calibrated, 
to within a relative accuracy of 10"*, to a photospheric radius of 
6.9599 X 10'" cm, a surface luminosity of 3.846 x 10" erg s"' and 
a ratio Zs/Xs = 0.0245 between the surface abundances by mass 
of heavy elements and hydrogen, at a model age of 4.6 Gyr. Fur- 
ther details on the model and oscillation calculations were given by 
IChristensen-Dals gaard ( 2008a, bj. 

To explore the sensitivity of the oscillations to the structure of 
the overshoot region we have considered a range of models with 
various characteristics, including models exhibiting the 'classical' 
sharp overshoot profile, a sequence of models converging towards 
the reference model, as well as models exhibiting various degrees 
of smoothness in the transition of the temperature gradient. Proper- 
ties of the models are provided in Table|2] In discussing the models 
with overshoot or otherwise modified, we use Model S as reference. 

Although the properties of the overshoot region are described 
in terms of the temperature stratification, the effect on the oscil- 
lation frequencies is predominantly controlled by the sound speed 
which in addition is affected by the composition. Of particular im- 
portance is the relatively steep gradient in the hydrogen abundance 
{X) established just below the fully mixed region which, as men- 
tioned above, includes the overshoot region down to r = rf. This 
affects the sound-speed gradient 

_ dine- „ din u 

V,2 = - ^V-—^, (16) 

d In p d In p 

where the last approximation used Eq. ([!}. This behaviour is illus- 
trated in Fig.|6] It is evident that the gradient in X a ccentuates the 
gradi ent in c" in Model S without overshoot (e.g., ISasu & Antial 
1 1994) . and similarly produces a step in Vj,2 at the edge of the over- 
shoot region in Model C2 illustrated in the figure. 

One might argue that further motion beyond the overshoot re- 
gion would cause additional diffusive mixing, thus likely smooth- 
ing out the step in the sound-speed gradient. To investigate the ef- 
fect of this on the signal we have computed Model OD, correspond- 
ing to Model C2 but with additional diffusive mixing in a small re- 
gion beneath rf; specifically, the maximum diffusion coefficient, at 
the edge of the overshoot region, was 50cm-/s, decreasing rapidly 
with depth. This causes a modest smoothing of the hydrogen profile 
which evidently is sufficient to eliminate the step in V^2 . 

As a further background for the discussion of the model sig- 
nals below, we first present a range of examples of profiles of 
the sound-speed gradient obtained with the formalism described 
in Section [2] These were chosen to provide an impression of the 
sensitivity of the signal to the detailed structure near the base of the 
convection zone, as well as to illustrate the structure that may be 
consistent with the observations. The input parameters of the mod- 
els, and some relevant characteristics, are listed in Table|2] Here, 
is the acoustic depth of the boundary r^z of the convectively unsta- 
ble region, while r^-iyi is the acoustic depth of the point of steepest 
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Figure 6. Properties of reference Model S (black solid lines), overshoot 
Model C2 (red dashed lines) and Model OD (blue dot-dashed line) which 
in addition has turbulent diffusion below the overshoot region. Top panel: 
Hydrogen abundance X against fractional radius rjR. Bottom panel: The 
thin curves show the temperature gradient V = din T/dlnp, and the heavier 
curves the coiTesponding gradient in the squared sound speed (cf. Eq. |16K 
As in Fig. [T]the vertical dashed lines mark the location of r{, ri and in 
Model CI. 
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Figure 7. Gradients V^2 of the squared sound speed in the reference 
Model S (solid line) and Models Al (dashed line) and A2 (dot-dashed line) 
with a sharp transition at the bottom of the overshoot region. 



© 2010 RAS, MNRAS OOO.fnfTS] 



8 7. Christensen-Dalsgaard et al. 



Table 2. Model properties. Here S refers to Model S of IChristensen-Dalsgaard et alj jl99^ . S' to a similar model, but with modified yi in the near-surface 
layers (cf Eq. lAlOl . and models Al - C3 were obtained utilizing the overshoot formulation described in Section |2] and characterized by ri, rf and yS, r^z and 
Tcz being the radius and acoustic depth at the base of the convection zone, defined by Vpad = V^d- Model OD in addition had a modest amount of turbulent 
mixing beneath the overshoot region while in Model OP the opacity near the boundary of the unstable region was reduced slightly (see text). To characterize 
further the properties of the models, T^^a shows the acoustic depth at the maximum in the overshoot region of dV/dr, and Cav/Hp is the total extent r^z - rf of 
the overshoot region, in units of the pressure scale height Hp at the base of the convection zone. Finally, i^o is the phase at zero frequency, determined from a 
fit to the eigenfunctions (cf. Appendix A). 
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Figure 8. Gradients V^j of the squared sound speed in the reference 
Model S (solid line) and Models Bl (dashed line), B2 (dot-dashed line) 
and B3 (triple-dot-dashed Une) converging to Model S (see also Table|2j. 



Figure 9. Gradients V^j of the squared sound speed in Models C 1 (solid 
line), C2 (dashed line) and C3 (dot-dashed line), illustrating increasingly 
smooth overshoot. 



gradient in V, providing an indication of the location of the acous- 
tic glitch. Also, 00 is the extrapolated near-surface phase at zero 
frequency, determined from a fit to the eigenfunctions (see Section 
lAlt . Finally, for comparison with the common characterization of 
overshoot. Coy I Hp shows the overshoot distance rcz - rf, in units of 
the pressure scale height at r^z. The properties of the models are 
illustrated in Figs|7]~[T0]in terms of V^2. 

Figure |7] shows that the present formulation emulates the 
sharp transition predicted by non-local mixin g-length theories 
(e.g.,IZahnlll99lh. a nd discussed in earlier work dBasu et al 111 994 
Basu & Antialll994lMonteiro etaPl 1994 [Roxburgh & Vorontsovl 



1994IChristensen :£alsgaard e t al.l 19951) which provided stringent 



constraints on overshoot of this form. Figure |8]shows a sequence of 
models converging to Model S (see also the parameters in TableO. 
A sequence of models of gradually increasing smoothness is illus- 
trated in Fig.|9] Finally, Fig.[TO]shows a model with a very smooth 



transition. To ensure a reasonable sound-speed profile, in view of 
the helioseismic inference (see also Section I6.lt , we shifted the 
boundary r^z of the unstable region outwards by about 0.002/? rel- 
ative to Model S, through a localized reduction of the opacity by 7 
per cent, centred near the bottom of the convective envelope. 

The frequencies of all models have been fitted with Eq. dl It 
in order to obtain the four parameters (aj , ai, td, yd)- Over 235 fre- 
quencies for each model have been used in the fit and no noise has 
been added to the frequencies. The values of the parameters found 
are listed in the first columns of Table |3] 

The values for A2.5 (see Ea.ll5t were obtained from the fitted 
parameters ai and ai which, though not listed in the table, are illus- 
trated in Fig. [TT] The amplitude values depend mainly on the local 
conditions at the base of the convection zone and are insensitive to 
changes in the model that do not affect this layer. Therefore the dis- 
tribution of the models in Fig.[TT]is an accurate representation of 
the sharpness of the transition at the base of the envelope. We note 
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Table 3. Seismic parameters resulting from the fit of the signal in Eq. (TT) to the frequencies of models. The first three columns are the results for noise-free 
data while the other six columns provide the results for 500 Monte Carlo simulations of each model when the uncertainties from a solar 72-d set of frequencies 
are considered. 
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Figure 10. Gradients V^2 of the squared sound speed in the reference 
Model S (solid line) and Model OP (dashed line), with a very smooth tran- 
sition and a small opacity reduction near the base of the convection zone. 
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Figure 12. The symbols show the values of the amplitude A2.5 versus 
acoustic depth of the signal for all models listed in Table |2] obtained 
by fitting the noise-free frequencies (values listed in Table|3). 
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Figure 11. The stars show the values of the amplitude parameters ui and 
£12 of the signal for all models listed in Table[3] The dotted line corresponds 
to ai = 02- 



in particular that, formally, 02 arises from a discontinuity in the first 
derivative in the sound speed, while ai arises from a discontinuity 
in the second derivative. Thus one would expect 02 to be bigger 
relative to ai for sharper transitions, as is indeed the case, although 
both terms contribute for all the cases considered. The correspond- 
ing inferred values of the amplitude A2.5 for the models are shown 
in Fig.M 

These results confirm wh at was already found by 
IChristensen-Dalsgaard et al] i ll 9951) : in order to have an amplitude 
below what is found for the standard model S, a sub-adiabatic 
region within the proper convection zone (above the Schwarszchild 
boundary) is required. 



5.2 Impact of noise on the seismic parameters 

When noise is added, the inferred values of the parameters are af- 
fected. The values in the rightmost columns of Table [3] were ob- 
tained from fitting model data after adding solar noise. For each 
model, Monte Carlo simulations were run with 500 realizations of 
errors characteristic of a 72-d solar observational dataset. 
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Figure 13. Examples for three models (A2 - stars; C2 - crosses; OP - filled 
circles) of the values 02 and ai obtained from Monte Carlo simulations us- 
ing solar uncertainties from a 72-d set of frequencies. The solid and dashed 
curves indicate the corresponding average combined amplitude /I2.5 (cf. 
Eq.[T5) and the ±1 o" error interval on A25. The values {a\,a2) obtained 
for the noise-free frequencies are also indicated as open circles for the three 
models (see Table[3). 



Figure 15. The stars show the average value from Monte Carlo simulations 
for the models using solar uncertainties, while the error bars are 1 cr errors 
associated with those uncertainties, as given in Table|3] The dotted fine cor- 
responds to a I =02- The small filled circles show the fits to the individual 
observed 72-d datasets while the open circles are for 1 year datasets. The 
dashed lines show the ± 1 o" range of A2.5 around the mean value for the 
solar data (full line). 




5I 

2170 2180 2190 2200 2210 2220 2230 2240 

(s) 

Figure 14. Values of and for two models, Bl and OP, inferred from 
the Monte Carlo simulation, based on observed errors. Note the strong cor- 
relation between these two parameters, already shown for the solar data in 
Fig.[5] The open circles show the results of fitting the noise-free frequencies. 

Fi gure [T3l illustrates the effect on the individual amplitudes ai 
and (32 for three models. It is clear from the figure that the added 
noise causes a systematic shift in the inferred values of a\ and 02, 
with 02 tending to be shifted to lower values, essentially corre- 
sponding to an apparently smoother transition. Similarly, although 
the value of A2.5 is more robust, it suffers a reduction which, for 
the sharper transitions, may exceed 2 cr: this is shown in Table [3] 
and is apparent for Model A2 in Fig. [13] Even so, the results con- 
firm that, when properly calibrated for the presence of noise, the 
value of A2.5 can be used to probe the sharpness of the base of the 
convection zone. 

As is the case for solar data, the largest correlation between 
inferred parameters is that between and y^: the effect of noise on 
these two parameters for model data is shown in Fig. [14] The stan- 
dard deviations of the infeiTed values of the parameters, and their 



mean values, are listed in Table[3] In this table, the values found for 
the seismic parameters in the presence of observed uncertainties are 
listed in the last 6 columns, together with the 1 a uncertainties ob- 
tained through Monte Carlo simulations of 500 realizations. These 
use the uncertainties for a 72-d set of solar frequencies reported 
before. 

As one would expect, the precision of the fitting is degraded 
when the amplitude is lower. In particular for and 7^, as lower 
amplitude decreases the ability to separate the frequency and mode 
degree contributions for the argument of the signal, leading to a 
wider variation, within the correlation found, of these two parame- 
ters (see also Fig.lI4t. 

We also note t hat, contrary to our pre v ious work 
jMonteiro et alj 1 1994 IChristensen-Dalsgaard et all Il995h . the 
option to estimate a priori the value of and fix the value of 
has resulted in an increase of the expected precision of the fitting 
for the three quantities listed in Table [T|(ai, 02, ta), required for 
testing the stratification of the overshoot layer This improvement 
is crucial in securing a proper calibration of the key parameters 
associated with the sharpness and location of the base of the 
convection zone. The increase in precision comes at the expense 
of systematic shifts between the noise-free fits and the averages 
of the model fits including noise, as is evident in Table [3] as well 
as Figs [T3] and [14] To correct for this the following comparison 
between the observations and the model results is based on the 
outcome of the Monte Carlo simulations for the latter. 



6 RESULTS AND DISCUSSION 

Based on the comparison of the fitting of the signal for the models 
(particularly in the presence of noise) and for the Sun, it is now 
possible to establish what are the implications for the stratification 
of the transition layer at the base of the convection zone. 
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Figure 16. The small filled circles show the fits to the individual 72-d 
observed datasets while open circles are for 1 year datasets. The stars show 
the average values obtained for Monte Carlo simulations of the models with 
1 o" en'or bars, as listed in Table [5] 



6.1 Comparing the models with the Sun 

Figure [TS] compares the amplitudes (aj and ai) obtained from so- 
lar data with those obtained from model frequencies with solar-like 
noise added (cf. the noise-free case in Fig. lilt . The model that 
seems most consistent with the solar observations, in this repre- 
sentation, is Model C2. This is confirmed by a comparison of the 
inferred values of A2 5 and for solar data and noisy model data 
(Fig. 1161 cf. the noise-free case in Fig. I12t. Model S and all the 
overshoot models in the A and B sequences exhibit amplitudes that 
are larger than what is inferred for the Sun. The models in the C se- 
quence seem to span the range of solar amplitudes, with Model C2 
providing the best fit. Model OP has a smaller amplitude than is 
inferred for the Sun, implying that the stratification in this model is 
even smoother than the average spherically symmetric radial strat- 
ification in the Sun. 

It is more problematic to draw inferences from a comparison 
of model and solar values of the acoustic depth parameter fj, be- 
cause the acoustic depth of the base of the convection zone is also 
strongly influenced by the sound-speed stratification in the rather 
uncertain near-surface layers. Specifically, the total acoustic size of 
the model and the acoustic location of the upper reflection bound- 
ary for each mode will have an impact on the value found for this 
parameter. Therefore its conversion into radial distance from the 
centre is not straightforward. We recall the values of obtained 
for Models S and S , shown in Fig. [T2l with a shift in the fitted 
values of fd (error free data) of about 25 s. The only difi'erence 
between these two models occurs at the surface, in order to bring 
the frequencies of S much closer to the solar frequencies. Conse- 
quently this value is an indication on the accuracy we may expect 
onfj. 

The fitting of the signal in the frequencies stills sufi"ers from 
some shortcomings, namely the incomplete representation of the 
dependence of the amplitude and argument of the signal on mode 
frequency and degree. On the other hand, the approach we have 
considered here uses of a large number of modes (over 240), by 
including data up to mode degree of 20. Also we isolate the signal 
directly in the frequencies avoiding noise amplification that occurs 
when the signal is isolated in frequency differences. Both aspects 
are unique to our approach allowing in the case of the Sun a di- 




Figure 17. Relative differences Sc /c~, inferred by means of SOLA in- 
version, between the Sun and three solar models, in the sense (Sun) - 
( model). The a nalys is used the so-called 'Best Set' of observed frequencies 
of iBasu et alj )l997 ). The vertical bars show the 1 cr eiTors in the dift'er- 
ences, estimated from the errors in the observed frequencies, while the hori- 
zontal bars provide a measure of the resolution of the inversion. In both pan- 
els the open symbols and solid line shows results for Model S. The closed 
symbols and dashed Une are for Model A2 (upper panel) and Model C2 
(lower panel). 



rect and very precise study of the base of the solar convection zone. 
However, it is evidently also of great interest to investigate the ex- 
tent to which information about the base of convective envelopes 
can be obtained from just low-degree data, as avail able in the fore- 
seeable future from observations of other stars (see lMonteiro et al] 
i200Q) . 

In this investigation we are using seismic data to look at only 
one aspect of the solar structure, namely the sharpness and shape 
of the transition in stratification near the base of the convection 
zone. The seismic data contain much additional information about 
the solar interior, notably about the run of adiabatic sound speed 
with depth. The sound-speed difi'erences between the Sun and three 
of the models considered here - Models S, A2 and C2 - inferred 
from helioseismic inversion of the data are illustrated in Fig. [17] 
The analysis was carried out u sing the so-called SOLA technique 
(see lRabello -Soares et alj 19991 for details on the imp lementation), 
foUowinglchristensen -Dalsgaard & Di Maurjl 20071) and using the 



1993) extending to 



combined 'Best Set' of frequencies (Basu e t al. 
/ = 99. Although there are statistically significant difi'erences in 
sound speed between the Sun and Model S, these difi'erences are 
small (0.4% in sound-speed squared, i.e. 0.2% in sound speed). The 
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differences between the illustrated models (A2 and C2) and the Sun 
are similarly small. In fact, although the models in this paper span a 
large range of behaviours in terms of the stratification near the base 
of the convection zone, we stress that they are all close to the Sun 
in terms of their small absolute differences in sound speed. 

We remark that the details of the sound-speed gradient arising 
from overshoot in some cases may have substantial effects on the 
seismically inferred sound-speed differences between the Sun and 
the model. In Model C2, illustrated in the lower panel of Fig.llVland 
in Fig. [6] the region where is larger than in Model S leads to an 
increase in the sound speed just below the convection zone which 
largely eliminates the bump in 6c^ found with Model S. However, 
if this overshoot region were to be shifted to greater depth, thus 
eliminating the region of sub-adiabatic gradient in the lower part 
of the convection zone, the effect on the sound speed would be 
amplified to such an extent as to lead to a strong negative sound- 
speed difference between the Sun and the model. (To avoid this 
for the very smooth Model OP we had to invoke a small opacity 
reduction to shift the location of the boundary of the convection 
zone.) Thus it remains important to test the models both based on 
the properties of the signal associated with the acoustic glitch and 
based on the inferred sound-speed differences. 



6.2 Implications for the stratification at ttie transition 

As found by IChristensen-Dalsgaard et al. I( ll995h standard models 
of the Sun either without or with the classical representation of 
overshoot, based on a non-local mixing-length formulation, are 
marginally inconsistent with the seismic data. In Fig. [16] both the 
amplitude for some models and for the solar data is shown. The 
reference Model S is outside the error bar of the data as well as any 
of the overshoot models calculated in the classical way. Even with 
a radiatively stratified overshoot layer the amplitude is found to be 
above the value found for the Sun. This is a clear indication that the 
transition at the base of the envelope requires a representation that 
is smooth over length scales of the order of the wavelength of the 
modes. 

This clear preference is a strong indication in favour of strat- 
ifications that include a sub-adiabatic layer within the proper con- 
vection zone, as a necessary condition to obtain sufficiently smooth 
transitions that will be consistent with the solar data. Any model 
that includes overshoot will need to be extended from within the 
convection zone, replacing the commonly adopted procedure of 
adding an overshoot layer to the convective envelope as defined 
by the mixing-length theory (or equivalent). 

The strong helioseismic preference for rather smooth profiles 
puts very strong constraints on many of the models we discussed in 
Sect.[T] In the following discussion we shall try to link the smooth- 
ness of the transition in the overshoot region to the underlying phys- 
ical processes. 

A common physical process relevant in most overshoot mod- 
els is buoyancy braking, which decelerates downflow plumes once 
a positive temperature perturbation has been built up after travelling 
downward for a sufficient distance in a subadiabatic stratification. 
Buoyancy breaking establishes a direct connection between the ki- 
netic energy of downflows at the base of the convection zone and 
thermal properties of the overshoot region beneath. This leads to 
a relation between the Mach number Ma required at the base of 
the convection zone, the subadiabaticity as well as the depth of the 
overshoot region. Assuming for simplicity a constant value for the 
superadiabaticity (5 = V - Va,], we can estimate temperature and 



density perturbations as function of the overshoot depth, z, as fol- 
lows: 



T' 

Y 



Q 



(17) 



Buoyancy braking gives a relation between the kinetic energy of the 
plume at the base of the convection zone and the buoyancy work in 
the overshoot region: 



-QV ~ -ggz= -gSgHp I — 



(18) 



where g is the gravitational acceleration. Combining both equations 
leads to (neglecting factors of order unity): 

M^.(^)\ (19) 

In the traditional non-local mixing-length overshoot models we 
have z - O.IH,, and 6 - 10'^ leading to M„ ^ 2 x 10"*, which cor- 
responds to a convective velocity of about 40 m s"' , consistent with 
mixing-length models assuming a filling factor of about 0.1. Hav- 
ing a smooth transition of the overshoot from nearly adiabatic to ra- 
diative temperature gradients requires that the average values for 6 
are of the order of 0.05. In that case our estimate requires - 0.05 
corresponding to a convective velocity of about lOkms"' at the 
base of the convection zone. The latter is (if at all) only possible if 
the filling factor for downflow plumes is very small. This is a in- 
evitable consequence of buoyancy braking, which applies to over- 
shoot models in a very general sense, whether numerical 3D simu- 
lation or any other simplified treatment. The key for the argumenta- 
tion above is the fact that a smooth overshoot profile implies auto- 
matically substantial deviations from adiabaticity in an extended 
region that is convectively mixed. Many of the numerical over- 
shoot simulations to date show smooth transitions. However, the 
latter are to a large degree related to numerical setups that have ei- 
ther a reduced stiffness in the radiative interior or have large energy 
fluxes (to reduce thermal relaxation time scales) that result in much 
larger Mach numbers and allow consequently for smo oth overshoot 
rofile s with substantial deviations from adiabaticity. iKapvla et al.l 
20071) presented a series of numerical experiments that show the 
return to more adiabatic step-like overshoot profiles with reduction 
of the overall energy flux. Extrapolating these results to the solar 
energy flux seems to make a step-li ke quasi-adiaba tic overshoot 
almost inevitable. The investigation o f lRempell ( l2004h aimed at un- 
derstanding the physical parameters under which a smooth transi- 
tion can be expected. The two necessary conditions found were that 
a) the downflow filling factor is small and b) there is a continuous 
distribution of downflow strength. While the latter is a natural con- 
sequence of convection, the former is very controversial and not 
supported by any numerical simulations we know about (the re- 
quired downflow filling factors would be smaller than 10"*). 

Overall it appears that within the framework of models dis- 
cussed above it is almost impossible to obtain sufficiently smooth 
profiles that agree with the helioseismic constraints. This difficulty 
is strongly related to the effective role of buoyancy braking in the 
overshoot region. Only if we relax this assumption is it possible 
to obtain smooth overshoo t profiles for mode r ate M ach numbers. 
It has been suggested by iPetrovav & MarikI ( 1 19951) that smooth 
overshoot profiles can be easily obtained once the assumption of 
a strong correlation between vertical velocity and temperature fluc- 
tuation is dropped. They suggested that in the extreme case of van- 
ishing correlation the penetration is only limited by turbulent dissi- 
pation, with no impact on the thermal stratification. For a smoothly 



© 2010 RAS, MNRAS QOQ.fTlfTsl 



Overshoot in the Sun 



13 



diminishing correlation a smooth transition of all overshoot quan- 
tities results^ 

IXiong & Dend jlOOlh , in fact, presented a model that leads 
to such very smooth transitions, by using a non-local convection 
model that is based on a moment approach computing a uto and 
cross -coiTelations of turbulent velocity and temperature l lXiongI 
Il989l) . This model predicts a strong decorrelation of velocity and 
temperature fluctuation in the lower overshoot region of the Sun, 
while the authors find a strong anticorrelation in the upper over- 
shoot region above the photosphere. The difference between the 
overshoot regions is attributed to the values of the Peclet num- 
ber, which are much smaller than unity in the photosphere and 
much larger in overshoot region at the base of the convection zone. 
The overall outcome is a temperature profile that is very simi- 
lar to the best match we found in this investigation. Their over- 
shoot region extends about 0.6//,, beneath the convection zone and 
shows already a significantly subadiabatic strati fication within the 
lower parts of the convection zone. Recently iMarik & Petrovavl 
(I2OO2I) used a different non-local model based on the formulation of 
ICanuto & Dubovikovl ( il997h and lCanuto & Dubovikovl j 19981) and 
found overshoot of only 0.06//p. While still much smoother than 
the sharp profiles of non-local mixing-length theory, this model 
provided again a temperature profile that stayed closed to adia- 
batic for significant extent of the overshoot depth. This class of 
models involves a number of assumptions, including the treat- 
ment of closure, and parameters generally determined from lab- 
oratory experiments or atmospheric or oceanographic measure- 
ments; thus a thorough study of their applicability under stellar 
conditions is indicated. In addition, a further investigation i s re- 
qui red of the differences betwe en the IXiong & Dend ( 1200 ll) and 
the lCanuto & Dubovikovl i ll 9971) approaches. 

Another possibility could be a spatial or temporal inhomo- 
geneity of overshoot leading to an average profile looking more 
smooth; however as estimated below this possibility seems unrea- 
sonable, too: 

Assume that we have a large filling factor ~ 0.1 so that the 
overshoot profile has a very sharp transition. Intermittent down- 
flows are able to perturb the sharp boundary between the almost 
adiabatic overshoot and the underlying strongly subadiabatic ra- 
diative zone. The consequence are buoyancy oscillations with the 
Brunt Vasaila frequency ojbv = ^jg\6\IHp. Due to the sharp transi- 
tion toward the radiation zone, oj^m would be given by the radiative 
values just beneath the overshoot region. With a value of dVjad/dr = 
10"^ m"' at the base of the convection zone, overshoot with a depth 
dos leads to a value oi 6 = -(dV|.ad/dr)(ios. Using d^s = 0.1 Hp, 
g = 500ms"^ Hp = 60 Mm yields w^v ~ 7 x 10"'* s"'. Amplitude 
d and velocity v of the buoyancy oscillations at the overshoot radia- 
tive interior interface are related by v = aiBwd, where v should be of 
the order of the convective motions perturbing the interface. Using 
a value of v = 100 m s"' yields d = 140 km, which is comparable to 
the thickness of the transition expected from the overshoot model in 
the first place. A larger amplitude requires a larger v, which can be 
only obtained by having faster downflows, meaning smaller filling 
factor (so nothing is gained here). We note that the above discussion 
implicitly assumed again a strong correlation between temperature 
and velocity fluctuations. 

It has been also speculated that the depth of overshoot depends 
on latitude, due to the c hanging influence of ro tation on convective 
motions. The study of iBrummell et al] 12002") found the deepest 
penetration at pole and equator and a minimum at around 30° lati- 
tude with about 60% of the maximum overshoot depth. This result 
was obtained by placing small Cartesian simulation domains (with 



periodic boundaries in the horizontal direction) at different latitude 
positions. Having a variation of overshoot depth in a global model 
adds additional complexity to the problem, since regions with dif- 
ferent vertical temperature structures have to be in horizontal force 
balance. These temperature differences can be estimated from Eq. 
( 1171 ) and can be quite substantial. Restricting T' to 100 K allows 
with 5 ~ 0.1, r ~ 2 X 10^ K and Hp = 50 Mm for only a displace- 
ment of 25 km. A global force balance in latitude would require 
either zonal flows or magnetic forces of substantial strength - at 
least the zonal flows should leave observable helioseismic signa- 
tures given the fact that only about 10 K temperature difference are 
required in the convection zone to balance the deviations of dif- 
ferential rotation from the Taylor-Proudman state. Without these 
flows or magnetic field such a configuration cannot be in a hydro- 
static balance and would return to spherical symmetry on a time 
scale given by cobv- 



1 CONCLUSIONS 

We have considered a variety of models with different stratification 
near the base of the convection zone, using a parameterization of 
the stratification that is motivated by models that have a spectrum 
of overshooting plumes with a range of depths of penetration. From 
the range of models we have considered, we find that a relatively 
smooth model (Model C2) fits the helioseismic data better than a 
standard model without convective overshoot, and better also than 
simple convective overshoot models which possess a more-or-less 
sharp transition to subadiabatic stratification beneath an overshoot 
region which is a nearly adiabatic extension of the convection zone. 
Thus a model characterized by an overshoot layer extending for 
iov/Hp - 0.37 is compatible with the solar data, corresponding to 
an additional fully mixed zone of about 3% in radius. Of course 
it is desirable that solar models be built whose stratification is de- 
termined self-consistently from the physics underlying the model, 
but the parameterization used here may be a reasonable way of ap- 
proximating such stratification in the context of simple spherically 
symmetric models. 

Our investigation method is \ery similar to that developed by 
iMonteiro et alj ( ll994h . though with some differences due to the fact 
that here we determine the phase ipg independent of our fitting to the 
oscillatory signal in the frequencies, and therefore we fit amplitude 
parameters ai and 02 separately. It should though be noted that we 
find that the fitted values of Oi and 02 are not independent, prob- 
ably because the frequency range of the modes used in the fitting 
is rather small and insufficient fully to separate the frequency de- 
pendent terms of the amplitude, and so their separate values cannot 
be relied upon. Indeed we find that noise in model data tends to 
cause the value of 02 to be significantly underestimated, relative to 
its fitted value for noise-free data. We find, however, that the com- 
bination amplitude A2.5 is a more robust indication of the sharpness 
of the transition in stratification near the convection-zone base. 

Our selection of models demonstrates that the signal analysed 
here is sensitive to acoustic glitches that would barely be noticed 
through helioseismic inversion which provides averages over such 
features. On the other hand, our analysis is by design insensitive to 
smooth differences between the Sun and the model, including very 
substantial differences in the sound speed which would be obvious 
from an inverse analysis. Thus these two techniques are obviously 
complementary in their ability to characterize the properties of the 
solar interior. 

It may be remarked that we have used models similar 
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to Model S, which incorpo r ates t he "old' solar element abun- 
dances of iGrevesse & Noelsl ( Il993h . rather than the "new" abun- 
dances suggested by the work of Asplund and collaborators (e.g., 
lAsplund et al.l 12003) . The new abundances are known to modify 
the opacities in the solar interior in such a way as to produce mod- 
els that are strongly in disagreement with the stratification inferred 
from helioseismology in the radiative interior, the largest discrep- 
ancy being in the vicinity of the base of the convection zone. Since 
our aim in this work is to investigate the subtle effects of overshoot 
in this region, it is prudent to start with models that are broadly 
consistent with the known stratification of the Sun's interior, rather 
than models based on the newer abundances that are further away 
from the Sun's stratification. 

Although the present analysis makes use of the full range of 
modes that are sensitive to the base of the convection zone, sim- 
ilar analyses are possible given just the low-degree modes avail- 
able in observations of distant stars in the foreseeable fut ure (e.g., 
iMonteiro et ai]|2000l . [Ballot et alj|2004l : IPiau et al.ll2005l) . This is 
potentially an important complement to the more detailed solar 
data, particularly given the very long timeseries and hence high 
frequency precision that will be obtained with the Kepler mis- 
sion (e.g., lGillilandetalJl2010l) as well as from planned dedicated 
ground-based facilities. 

From this investigation we can conclude that i) overshoot is 
necessary to improve the agreement between models and helio- 
seismic constraints, ii) the required overshoot profiles are outside 
the realm of the classic "ballistic" overshoot models and iii) the 
lower part of the convection zone is likely substantially subadia- 
batic. We cannot measure the latter directly, but our investigation 
indicates that the required level of smoothness cannot be achieved 
without iii). CuiTently only non-local convection models that are 
based on auto- and cross-coiTelations of velocity and temperature 
perturbations are capable of providing overshoot profiles with the 
desired degree of smoothness in the transition. However, these the- 
ories have hidden parameters and rely on closure models. To our 
current knowledge there has not yet been a thorough study on how 
robust the properties of overshoot (under the conditions found at 
the base of the solar convection zone) are within the framework of 
these models. We hope that our investigation motivates more re- 
search in that direction. The indication that the lower part of the 
convection zone could be substantially subadiabatic (with values 
of |V - V;,(]| > 10"') has profound consequences for the storage and 
stability of magnetic field and the overall role that the lower con- 
vection zone (not just the overshoot region) might play in the solar 
magnetic cycle. 
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APPENDIX A: DETERMINING THE PHASE OF SOLAR 
EIGENFUNCTIONS 

Al The eigenfunction phase 

The phase /) in Eq. ^ is obviously closely related to the 
eigenfunctions of the oscillations and hence can be determined 
from fits to those eigenfunctions. This was discussed in de- 



tail bvlChristensen-Dals gaard & P erez Hernande z ( 1992) (s ee also 
iRoxburgh & VorontsovllT996b . in the context of the buvalll ( fl982l) 
law. According to this, the frequencies satisfy 
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Figure Al. Scaled eigenfunction E = (pc)'^^r^r, for the mode with / = 0, 
n = 21 in Model S, as a function of acoustic depth r. The crosses show the 
results of fitting the asymptotic eigenfunction. 
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where L? = 1(1 +1). This is obtained from Eq. ([8} by applying the 
appropriate inner boundary conditions, with <p and a being related 
by 



a -\ — \ n . 
4 



(A2) 



For simplicity we neglect the higher-order terms in the near-surface 
behaviour and assume that a, and hence (p, are functions of a> alone. 
The quantity 0o entering into the fitting formula is the intercept at 
ai = in a linear fit to as a function of <i> and hence can be 

dete rmined from a similar fit to a. 

IChristensen-Dalsgaard & Perez Hernandez! d 19921) analysed 
the dependence of a on the properties of solar models by fitting 
Eqs {8} and ([9} to computed solutions of the equations of radial os- 
cillations in the models, assuming only the surface boundary con- 
ditions and hence determining a as a continuous function of a> ( a 
similar procedure has been used by [Roxburgh & Vorontsov|[l996h . 
Here we consider instead full eigenfunctions of the model, over a 
range of degrees, to approximate more closely the actual fits. As a 
simplification, we use only relatively low-degree modes and apply 
the fit rather close to the stellar surface. Then Eqs (|8) and (|9]( can 
be approximated by 



-1/2^-1 



co&{a)T + <p) . 



(A3) 



in terms of the acoustic depth t (cf. Eq. O. We therefore fit 
Acos{ll)t + cf>) to E = (pcY'^r^,. in a least-squares sense, over a 
suitable interval [ti, T2] in t. 

As an example we consider Model S 
jChristensen- Dalsgaard et all Il996l) . Figure lAll shows an ex- 
ample of the fit to an eigenfunction, for the mode with / = 0, 
n = 21. The fit was made in the interval between 500 s and 2000 s 
in T. As illustrated by the crosses, the fit is excellent. 

To determine (picj) we consider modes in the interval 5 < Z < 
10; also, in accordance with the mode sets used in the fit for the base 
of the convection zone, we consider frequencies between 1.9 and 4 
mHz. The phases resulting from the fit are shown in Fig. IA2I con- 
firming that the scatter at given frequency is modest. Also shown is 
a least-squares fit of a straight line which yields <po = 0.890, and a 
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Figure A2. Fitted phases for modes between 5 and 10 in degree and 1.9 
and 4 mHz in cyclic frequency v = Li>l2n, in Model S. The solid and (barely 
visible) dashed curves show linear and cubic least-squares fits to the points, 
respectively. 



Figure A4. The predominantly frequency-dependent residual difference T 
(cf. Eq. lA9t resulting from a fit to frequency dift'erences between Model C2 
and Model S, using Eq. )A8t . The straight line shows a linear least-squares 
fit to the results, resulting in a determination of 0o- 
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Figure A3. Duvall plot for the modes between 5 and 100 in degree and 1.9 
and 4 mHz in frequency, in Model S, determining a from the cubic fit to the 
phase in Fig. lA2l 

cubic fit which follows the computed points very closely but yields 
an intercept at toi = of 0.503. 

The Duvall relation resulting from (p, and hence a, determined 
in this manner is illustrated in Fig. lA3l using a{(jj) as obtained from 
the cubic fit in Fig. IA2I Here all modes with 5 < / < 100 and 
frequency between 1.9 and 4 mHz are included. The match to the 
Duvall law is very good, strongly indicating that as determined 
from the low-degree modes can be applied over a broad range of 
degrees. 

A2 Differential phase determination 

While the eigenfunction fit provides the most appropriate phase 
0(a)) for the models, and hence presumably an appropriate estimate 
of this is obviously not possible for the observations. However, 
the results above suggest that an estimate of and 0o can be ob- 
tained from a fit of the Duvall relation. One possibility would be 
to determine the function a{(jj) that most successfully collapses the 
points in the Duvall plot to a curve and to obtain ^{oS) from that. 



A possibly simpler approach is to use the diff erential form of the 
Duvall law dChristensen-Dalsgaard et al ]ll988h to determine a cor- 
rection to 00, relative to a suitable reference model. By linearizing 
Eq. JAU in small changes dc and 6a we obtain 



5^ = 

Cl) 


--'HiiojIQ + 'Hiioj) , 




(A4) 


where 
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'"1 LVr"'dr 
\ cuV] c ' 


da 
dci) 


(A5) 


•Hi(w) 




dr 
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(A6) 


and 










n 

= —5a(LL)) . 




(A7) 



As discussed bv lChristensen-Dalsgaard et al.l ( Il989l) the functions 
Hi and "7^2 can be approximated by a double-spline fit to the scaled 
frequency differences, to obtain 'Ki and 'Hi. Given "Hi, Eqs JA2t 
and l lA7t show that the difference in the phase can be obtained as 

(50 = -Ji6a ^ -w-Kz . (A8) 

It should be noticed that a fit of the form given in Eq. ( IA4t 
can only determine the functions Hi and Hi to within a constant. 
Consequently, S(j) as determined from Eq. (IA8) is only determined 
to within a constant multiple of a). Obviously, this has no effect on 
the inferred value of ScfiQ. 

We have tested this procedure by applying it to a few of the 
models discussed in Section lSTT] Here the frequency range was, as 
usual, between 1.9 and 4 mHz, and modes of degree between and 
100 were included. For simplicity, the term in dcy/doj was neglected 
in computing S from Eq. ( IA5t . To determine (50o we carry out the 
fit to determine the function "Hi (oi/L) and evaluate 

T = -{SSoj-uj'Hi)\ (A9) 

we then carry out a linear least-squares fit to T, as a function of 
frequency, resulting in (50o- The procedure is illustrated in Fig. lA4l 
for Model C2. In Table I All the result is compared with the actual 
difference in 0o, relative to Model S, obtained by fitting the eigen- 
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Table Al. Phases and phase differences resufting from eigenfunction and 
differential Duvall fits; i50q ' is the difference relative to Model S computed 

(2) 

directly from 0o for the two models, while Sf^ is the corresponding dif- 
ference infen'ed from the differential asymptotic Duvall fit (Eq. |A8K The 
solar result of the differential asymptotic analysis used a 72-d MDI dataset 
starting in April 1996. 



Type 


ID 


00 






Model 












S 


0.8904 








C2 


0.8944 


0.0040 


0.0037 




OP 


0.8954 


0.0050 


0.0047 




s' 


0.5485 


-0.3419 


-0.3257 


Sun 








-0.3181 



0.20 



w= 0.10 




0.00 
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Figure AS. The predominantly frequency-dependent residual difference T 
resulting from a fit to frequency differences between Model S and Model S, 
using Eq. jA8t . See caption to Fig. lA4l 



functions in the two models; the agreement is clearly quite good. 
Table lAll also lists the very similar results for Model OP. 

As a more extreme example, we consider a model. Model S , 
based on Model S but modified to suppress the dominant, near- 
surface part of the difference s between the observed frequencie s 
and those of Model S (e.g.. IChristensen-Dalsgaard et al.|[l996h . 
Specifically, Fj was modified by multiplying it by 



1 - 0.4 exp 



{r-roY 



Ar2 



(AlO) 



where ro = (1-2.5x10"*)/? and Ar = 1.5x10"*/?. This was de- 
signed to emulate the effect of turbulent pressure on the thermody- 
namics of the strongly superadiabatic part of the convection zone 
(C. S. Rosenthal, private communication). The differences between 
the frequencies of this model and those of Model S essentially 
match the corresponding differences for the observed frequencies 
and hence contain a large frequency-dependent component (after 
scaling). This is reflected in the large frequency-dependent part of 
the frequency differences and hence in the resulting T, illustrated 
in Fig. IA5I The value of (50o obtained from fitting to T is listed 
in Table lAll together with the value resulting from the fits to the 
eigenfunctions. Remarkably, even for this very substantial differ- 
ence there is reasonable agreement with the directly determined 
values of for the models, indicating that the differential Duvall 
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Figure A6. The predominantly frequency-dependent residual difference T 
(cf. Eq.|A9j resulting from a fit to frequency differences between MDI ob- 
servations and Model S , using Eq. ^A8). The observations were the 72-d 
set used in Table lAl] See caption to Fig. lA4l 



analysis provides a robust method for estimating 0o> also for ob- 
served frequencies. 

We have also applied the differential Duvall analysis to a 72-d 
set of MDI observed frequencies. As shown in Table lAll this yields 
a (500 very similar to that obtained for Model S ; this was to be ex- 
pected, given that this model was fitted to the observed frequencies. 

In the fits of model frequencies to Eq. dl It we used the value 
of 00 resulting from a fit to the corresponding eigenfunctions. In 
the analysis of the observed frequencies, we have used Model S 
as a reference; in this way we hope to minimize any systematic 
effects that would result from the larger frequency differences rela- 
tive to Model S. As an example, Fig. lA6l shows the phase-difference 
plot for one of the 72-d observational sets. The linear fit gives 
"500 = -8.9 X lO"*, again indicating that Model S provides a 
good fit to the observations. It is likely that the oscillatory com- 
ponent arises from a combination of the residual effects of the 
acous tic glitch associated with the second helium ionization zone 
(e.g. lGougrill990l ; Pv^rontsov et al]|l99ll; iMonteiro & ThompsonI 
l2005h and a component of the near-surface difference which has 
not been eliminated by the modification l lAlOt to Fi . The former 
effect would indicate that the helium abundance (or the equation 
of state) in the model (and hence in Model S) is not in full agree- 
ment with that of the Sun; taken at face value the variation below 
2500 /jHz suggests that the envelope helium abundance of 0.245 in 
the model is too low by about 0.008. 



A3 Solar-cycle variation 

We have carried out the differential fit to determine the value 
of (500 and hence 0o for all the datasets involved in the analy- 
sis. The results for the 72-d sets were illustrated in Fig. |2] which 
clearly reflects the variation in solar activity. To understand the be- 
haviour of the variation we compare observed frequencies at an 
arbitrary phase of the solar cycle with frequencies at solar mini- 
mum and recall that the frequencies increase with solar activity, 
the change being a steeply increasing function of frequency (e.g., 
iLibbrecht & . Woodard 1990). We also note that "Hi essentially cor- 
responds to the scaled relative frequency differences, apart from 
the arbitrary constant in the separation into "Hi and 'H2 in Eq. l lA4t . 
Thus increases with increasing frequency. Choosing the con- 
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Figure A7. The functions —U/H2 resulting from the spline fit (cf. Ea. lA41 
to the frequency dift'erences between all sets of 72-d observations and 
Model S , after subtraction of the linear fits to the corresponding T. 

slant such that 'H2 is zero at the low end of the frequency range 
considered, the function f is similarly zero at low frequency but 
with a negative slope. Consequently, the intercept (50o at zero fre- 
quency, which is obviously unaffected by the choice of constant, is 
positive, as observed in Fig. |2] and increases with increasing solar 
activity and hence increasing frequency. 

There have been suggestions that the effect of the glitch asso- 
ciated with the seco nd helium ioni z ation zone varies with the phase 
of the solar cycl e teough 1995, '2002b!; iBasu & Mandel |2004 
iBallot et al.ll2006 l: Ivemer et al. 2006). We would expect that such 
a variation should be visible also in T ■ To investigate this, Fig. lA7l 
combines the fitted 'H2 for all 72-d observations, after subtraction 
of the linear fit. There is evidently very little scatter in the result, 
and further inspection shows no evidence for systematic variations 
with the phase of the solar cycle. Similarly, the second differences 
in the smooth component of the frequencies (see the lower panel 
of Fig. [3) show no significant variation with solar cycle. To set the 
scale, we note that the change 5v in cyclic frequency, corresponding 
to a change ST in T , approximately satisfies 

Av 

5T . (All) 

n 

From this we estimate that the range of the residual variation in y 
with solar cycle, after removing the linear trend, is below 0.04 /jHz, 
which apparently is substantially smaller t han the variation in the 
sig nature of helium ion ization inferred by ISasu & Mandell Jiooj) 
and lVemer et al.l(l2006h . This deserves further investigation. 



